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Pesticide Precautionary Statement 


This publication reports research involving pesticides. It does not contain 
recommendations for their use, nor does it imply that the uses discussed here have been 
registered. All uses of pesticides must be registered by appropriate State and/or Federal 
agencies before they can be recommended. 


Caution: Pesticides can be injurious to humans, domestic animals, desirable plants, and 
fish or other wildlife—if they are not handled or applied properly. Use all pesticides 
selectively and carefully. Follow recommended practices for the disposal of surplus 
pesticides and pesticide containers. 


Foreword 


This report is published as a part of the USDA Forest 
Service program to improve aerial application of pesticides, 
specifically by using pesticides and delivery systems tailored 
to the forest environment. The program is conducted jointly 
by the Technology and Development Center, Missoua, MT , 
and the Forest Pest Management Staff, Washington Office 
at Davis,CA, under the sponsorship of State and Private 
Forestry. 


Details of the aerial application improvement program are 
explained in two Forest Service reports: A Problem 
Analysis: Forest and Range Aerial Pesticide Application 
Technology (Equipment Development Center Report 7934 
2804, July 2979, Missoula, MT) and Recommended 
Development Plan for An Aerial Spray Planning and 
Analysis System (Forest Pest Management Report FPM 82- 
2, February 2982, Davis, CA). 


The analysis was cosponsored by the Physics Division— 
Research Directorate—(Chemical Research, Development 
and Engineering Center)—U.S. Army . 


The study was conducted as part of Program WIND (winds 
in nonuniform domains). Authority for the cooperative 
program is the Supplemental Agreement, dated February 
1985, to the master Memorandum of Understanding 
Between U.S. Department of Defense and U.S. Depa” ant 
of Agriculture Relative to Cooperation with Respect to © 20d, 
Agriculture, and Other Research of Mutual Intesrest. 
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EXECUTIVE SUMMARY 


Ground deposition predictions from the computer program AGDISP are compared 
with four helicopter trials from Program WIND Phase III (EMCOT). In addition, the 
simplified downstream wake flow field assumed in AGDISP is compared with a detailed 
free wake prediction of the trailed .helicopter rotor flow field. Results indicate the 
following: 


1. AGDISP does a reasonable job of predicting the shape of the ground 
deposition pattern in Program WIND. 


2. The neutral atmospheric background model in AGDISP should be extended 
to nonzero Richardson number to modify the ambient turbulence level and 
improve predicted deposition in stable and unstable layers similar to those 
present in Program WIND. 


3. The Program WIND Phases I and III and Withlacoochee test results should 
be re-examined now that AGDISP contains a canopy model. 


4. Additional tests on flat open terrain with nominal wind speeds are strongly 
recommended to complete validation of AGDISP helicopter ground 


deposition predictions. 


5. Comparison of the AGDISP helicopter wake flow field model with a more 
complicated free wake computer model shows excellent agreement between 
the two models and supports the continued use of the simplified AGDISP 
model. 
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1. INTRODUCTION 


The computer program AGDISP predicts the motion of aerially released material 
including the mean position of the material and the position variance about the mean as a 
result of turbulent fluctuations (Teske, Ref. 1). AGDISP is based on a Lagrangian 
formulation of the released material equations of motion, and includes simplified models 
for aircraft wake and ambient turbulence effects. To this point AGDISP predictions have 
been compared successfully with numerous fixed-wing aircraft data sets (Teske, Ref. 2 for 
example). Comparisons of the helicopter model have been attempted in the following 
(leskeeneta): 


1. Forest Service test results from the Withlacoochee seed orchard in 1980 with a 
Hughes 500C helicopter; 


2. Program WIND Phase I test results from the Chico almond orchard in 1985 
with a Hiller 12E helicopter; and 


3. Program WIND Phase III test results from the Red Bluff forest in 1986 with a 
Bell 206 helicopter. 


These tests were all conducted in a forest canopy; these comparisons were all 
performed before AGDISP contained a canopy model. Data taken in an open field would 
remove any canopy influence and provide a clearer data base for AGDISP comparisons. 
The Program WIND Phase II (EMCOT) data set, collected in open sloping terrain (Ekblad 
et al., Ref. 4) satisfies many of the requirements for such a study. To this end the four 
reduced trials (numbers 8, 9, 10 and 11) are compared with AGDISP predictions in Section 
2 of this report. 


Since its inception, the assumed helicopter wake flow field model in AGDISP has 
never been compared with data. It appears likely that it never will, since data in close 
proximity to a helicopter rotor blade is either rare or proprietary. An alternate approach, 
comparing the assumed AGDISP wake model to the flow field generated by a far more 
complex (and far more physically detailed) computer model, is examined in Section 3. 
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2. GROUND DEPOSITION COMPARISONS 


Program WIND Phase II] (EMCOT) Spring Tests were conducted from late April 
to early May 1986 at a cleared sloping site in the Sierra Nevada foothills approximately 25 
miles east of Red Bluff, CA. This area consists of a pair of ridges, between which runs a 
small stream (Figure 2-1). An anemometer tower grid recorded wake velocities of aircraft 
traversing normal to the grid. Three card rows were set along the two ridges (Rows A and 
C) and in the gully at the base of the towers (Row B). Additional meteorological stations 
around the test site recorded local wind conditions. 


All three card rows had a common origin, at A-0, B-0 and C-0, as seen in Figure 
2-1. This position is designated 0.0 meters in all that follows. Positive distance is 
downslope to the west. The aircraft flew north to south along lines parallel to the flight line 
shown in Figure 2-1. The ground slopes away from the right wing of the aircraft (average 
downslope is about -3.5 deg). 


In all, 49 aircraft flyovers were recorded on the anemometer tower grid (Teske, 
Ref. 5). Of these runs, 26 flyovers included the collection of ground deposition data along 
the three card rows. All of this data were subsequently reduced by the USDA Forest 
Service (R. B. Ekblad, private communication). Analysis of this data and the tests 
themselves point out the following sources of error: 


1. Ekblad et al. (Ref. 4) imply that the assumed spread factor relating card 
stain size to actual drop size may be in error by up to 20 percent. This 
suggests that the evaluated mass on the cards could be off by as much as 50 
percent. 


2. Post-test inspection of the helicopter spray system detected leaks that 
generated extraneous large drops (above 600 pm), presumably at the 
expense of smaller drops (below 600 Lm). 


3. Observations at the test site indicate that the helicopter pilot may have started 
the spray too late and/or stopped the spray too soon in several of the trials. 
In these cases the local winds could have drifted the spray beyond a card 
row before ground deposition. 
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Figure 2-1: | Program WIND Phase III (EMCOT) Red Bluff test site. 
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Integration of the ground deposition pattern from the reduced card data reveals that 
in many trials, less than one-half of the sprayed material reached the ground. The 
ambiguity generated by the above three points suggests that any predictions by AGDISP 
should be normalized by the actual fraction of material recorded on the card rows. This 
Operation makes this test series less than perfect when comparing against AGDISP 
predictions. Nonetheless, it seems reasonable, since AGDISP will account for all released 
spray material. 


Four trials (numbers 8, 9, 10 and 11) were conducted between 6 and 8:30 am on 
April 30th. This time period forms the basis for a detailed analysis of meteorological 
conditions (Hauser et al., Ref. 6) and provides the only trials for which surface winds have 


been determined. The time-averaged crosswind velocity profiles are shown in Figure 2-2. 


AGDISP predictions were made for each of the three card rows for each of the four 
trials. Helicopter characteristics are summarized in Table pil Specific trial and card row 
information is presented in Table 2-2. Each AGDISP run combines six drop sizes as given 
in Table 2-3. The sizes are selected consistent with the card deposition data reduction 
algorithm (R. B. Ekblad, private communication). The mass fraction for each size was 
obtained by summing the drop size distribution data from a wind tunnel test of the nozzle 
and spray material (Yates, Ref. 7). AGDISP comparisons with data are shown in Figures 
2-3 to 2-14. A typical AGDISP input deck is shown in Table 2-4. 


The "Deposition Fraction" column in Table 2-2 illustrates the material fraction of the 
0 to 600 um drops measured on the cards. Trial 9 Row C is particularly surprising, and 
highlights the difficulty of using this data set for a definitive comparison with AGDISP 
predictions. Once AGDISP predictions are multiplied by the "Deposition Fraction," 
however, the comparisons with data show great similarity for both the Liters per Hectare 
and the Drops per Square Centimeter deposition plots. The best agreement is Trial 11 Row 
B in Figure 2-12. The problem with Trial 9 Row C is reflected in the Drop plot in Figure 
2-8; a similar problem exists in Trial 10 Row C of Figure 2-11. Most of the rest of the 
comparisons all show the prediction off to the left of the data. Two reasons may be 
identified for this discrepancy: 
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Figure 2-2: 


Crosswind velocity profiles for card rows A, B and C for trials number 8, 
9, 10 and 11. 





TABLE 2-1 
Program WIND Phase IJ (EMCOT) Test Characteristics 


Aircraft Type Bell 206 Helicopter 
Nominal Weight 9342 N 
Rotor Radius 5m 
Blade Rotation Rate 390 rpm 
Number of Nozzles 22 
Application Rate 19.4 gal/min 
Atomizer D8-45 straight back 
Tank Mix 67% water 

33% Glycerine 


2 lbs/gal No. 2 Blue Food Coloring 
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TABLE 2-2 


Program WIND Phase II] (EMCOT) Card Row Characteristics 


Trial Card 


Aircraft 


Offset Aircraft Deposition 


Number Row Height(m) Distance (m) Speed (m/sec) Fraction 


8 A 14.9 
B 22.4 
GC 17.0 
2) A 13.4 
B 20.9 
Cc 1553 
10 A 127 
B FAVS 
e 14.6 
i A 16.0 
B 2305 
Cc 18.1 
Aircraft Height (m): 
Offset Distance (m): 
Aircraft Speed (m/sec): 


Deposition Fraction: 


51.8 eal) 0.498 
SES 220) 0.329 
Saige 25.0 0.447 
15S 22.4 0.282 
51.8 22.4 0.289 
51.8 22.4 0.076 
B32 20.6 0.354 
oa e8) 20.6 0.287 
Sa 20.6 0.201 
Ol 22.4 0.276 
AUEt 22.4 0.230 
TAG 22,4 0.206 


Apparent height of the aircraft over the card row at the 
Offset Distance. 


The distance from the zero line (A-0, B-0 or C-0) to the 
aircraft flight line. 


The speed of the aircraft. 


The mass fraction of material collected on the cards for 


drop sizes between 0 and 600 Lm, divided by the 
anticipated mass fraction from wind tunnel test data 
(Yates, Ref. 7) and the application rate. 
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TABLE 2-3 


Program WIND Phase III (EMCOT) Size Distribution Characteristics 


Nini 
Diameter ([1m) 


200 
300 
400 
500 


Maximum Average Mass 
Diameter ({1m) Diameter (1m) Fraction 
100 63.0 0.0185 

200 155.4 0.1029 

300 2303 0.1701 

400 352.4 0.2368 

500 451.8 0.2132 

600 55155 0.1470 


Average Diameter (Um): 


Mass Fraction: 


Computed from volume average formula in 
Herdan (Ref. 8). 


Summed and interpolated from wind tunnel test 
data (Yates, Ref. 7). 
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Figure 2-3: 




















AGDISP composite ground deposition predictions (solid curves) compared 
with card data (squares) for Trial number 8, Row A. The top figure is mass 
in Liters per Hectare; the bottom figure is number density in Drops per 
Square Centimeter. 
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Figure 2-4: 








Trial number 8, Row B. 
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Figure 2-6: Trial number 9, Row A. 
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Figure 2-9: Trial number 10, Row A. 
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Figure 2-12: Trial number 11, Row A. 
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Figure 2-14: Trial number 11, Row C. 
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TABLE 2-4 
Sample AGDISP Input File 


Q000 EMCOT TRIAL( 11 LINE ’B 
O07 OFS U0 Oa 

OOD a LO Use Owes 0 

0015" 30 


0020) 94 wo UZ joe 2 4 
0029 0.0 Den O20 
0029 Zine Seem) Sales 
0029 10 1 a 
0029 dy pare Zhe eve: 
0029 235 25 =O118 
COZ eS fod 320 =O 70S 
0030 342:,0873 907.0 


0050 S0ReOrO 

OLONSOT Rorat TO Oe aA een OO MMe tet Oe ea 8) 
OVO Le =52.03 2 
00Gb C457 =F 20816 
DOG I a4. 2 Va VES) 
UVUG I =3:-06 =25.10 
OO6iiet EBi2 T er. h6 
OO 6 1g 214 LO 
OG ire =257 10 
0061 Mai .63 rer eno 
OO G1 ae = Ae lO) 
CU Or genes es) doen Ze 
OCG GOW Sve K-21. 76 


0061 O55) me 266 
0061 04 Gee Zan 0 
0061 mice 2288 6 
0061 L403 Ah 
0061 P26 ae Sy ae TRS 
0061 2.74 2h eh 
0061 3x21 aa eG 
0061 B00 =n 
0061 AG 2 FF AO 1.6 
0061 4a, eal ae IBS 
0061 aes RNS 


fag) 





1. The time-averaged crosswind velocity profiles (Figure 2-2) may not 
faithfully represent the actual winds during deposition. Previous AGDISP 
comparisons (Teske, Ref. 2) show excellent agreement with ground 
deposition data when crosswind is present. The crosswind velocity profiles 
are imposed on the AGDISP trajectory solution; if these profiles are 
incorrect, AGDISP cannot generate accurate ground deposition plots. 


2. The position of the aircraft relative to A-0, B-0, C-0 may be incorrect. Field 
observations provided this information. 


These trials were conducted in the early morning hours when a temperature 
inversion and accompanying stable boundary layer are generally present. The AGDISP 
model considers only neutral conditions; the presence of a stable atmosphere would 
decrease the effect of ambient turbulence on the growth rate of the spread of the released 
material about the mean material trajectories. Atmospheric stability is a function of the 
Richardson number, a nondimensional parameter that relates the local temperature gradient 
to the local velocity gradient (Monin and Yaglom, Ref. 9). Simplified theories of 
turbulence (Bilanin et al., Ref. 10) may be interpreted to compare the effect of stability on 
ambient turbulence level. Such a plot is given in Figure 2-15. For the stable conditions of 
the present study, the anticipated Richardson number would be positive, the "Turbulence 
Factor" would be less than one, the growth rate of the spread of the released material would 
be diminished, and the predictions in Figures 2-3 to 2-14 would all exhibit larger, narrower 
depositions as suggested by the data. The effect of stability on turbulence should be added 
to the AGDISP computer program, and could possibly be implemented with a stability ratio 
(SR) approach similar to that suggested by Akesson and Yates (Ref. 11) to characterize 


ambient atmospheric conditions. 
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3. WAKE FLOW FIELD COMPARISONS 


The helicopter wake flow field model in the AGDISP computer program (Teske, 
Ref. 1) includes both the hover downwash and the generated tip vortex pair by partitioning 
the helicopter weight between the two effects as a function of distance downstream. The 
hover downwash model is taken from actuator disk theory for a propeller, while the 
strength of the tip vortex pair is found by assuming a rectangularly loaded rotor blade 
circulation distribution. The partitioning is done by an exponential function F that begins 
at 0.0 (the wake is all helicopter downwash) and quickly increases in value toward 1.0 (the 
wake is all vortex pair). A model constant in F was set by comparison with preliminary 
free wake forward flight code predictions (Bliss et al., Ref. 12). 


The simplified AGDISP model for the growth of the tip vortex circulation strength 
may be compared with predictions from the most detailed (and available) helicopter wake 
flow field model (Quackenbush and Bliss, Ref. 13). This code utilizes a numerical 
analysis based on the curved vortex filaments described in Bliss et al. (Ref. 12). The 
model is composed of vortex filaments laid down along contours of constant sheet strength 
in the wake. The skewed/curved filaments provide a natural representation of the free wake 
vorticity field, which simultaneously accounts for both shed and trailed vorticity. An 
additional advantage of the model is that it provides a visually meaningful representation of 
the wake since the filaments correspond to the actual resultant vorticity field. Vortex 
filaments leaving the blades are all of constant strength and equal value. For this reason 
close spacing between the filaments implies a strong net influence from that region of the 


wake, whereas a sparse spacing indicates a region having little effect. 


Figure 3-1 presents the results of running this computer model for the Bell 206 
helicopter conditions given in Table 2-1 previously. Figure 3-1 shows the vortex filament 
structure trailing off one of the two blades (for clarity). It is quite complicated and, when 
combined with filaments from the other blade, generates the wake flow field downstream 
of the helicopter. Crossplane velocity vector plots are shown at one, two and three blade 
radii downstream in Figures 3-2a, 3-2b and 3-2c respectively, compared against AGDISP 
helicopter model predictions for the same initial conditions. The familiar vortex-pair-like 
velocity pattern is very evident. So too is the anticipated lateral asymmetric structure of the 
downwash flow field in the complicated helicopter model results. The nght (advancing 
blade) side exhibits a somewhat weaker tip vortex than the left (retreating blade) side. At 


the later radii the mean position of the vortex pair is seen to move downward. 
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Figure 3-2a: Crossplane velocity vectors at a downstream distance of one radius (top) 
compared with AGDISP predictions (bottom). 
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Figure 3-2c: Crossplane velocity vectors at three radii downstream. 
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The complicated helicopter model calculation took six hours on a MicroVAX II 
computer; a typical AGDISP run takes one to two minutes. 


To reproduce the size of the velocity vectors in the complicated helicopter model 
required the imposition of a cutoff distance in the AGDISP helicopter model of 0.1R. This 
minor revision to the AGDISP model will reduce slightly the influence of the vortex when a 
droplet is very close to a vortical center. Further out, there will be no change in the 


influence of the vortex on the motion of material droplets in the simulation. 


Recent anemometer tower data reduction for Program WIND Phases I and III 
(Bilanin et al., Ref. 14) lead to the refinement of a least squares influence coefficient 
technique for recovering the position and strength of the vortex pair generating the 
anemometer signals. This effort was directed toward discerning the decay rate of vortex 
pairs in the atmosphere, but the same technique can be just as easily applied to the 
crossplane velocity vectors shown in Figure 3-2 to recover their position and tip vortex 
circulation strength. The growth of this circulation strength may then be compared with the 
AGDISP model, giving the surprisingly good result shown in Figure 3-3. The exponential 
wake model in AGDISP does a remarkably good job of simulating the complicated wake 
flow field generated by the far more detailed helicopter code. 
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